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In this study, based on the self-energy method and the total energy calculation, the indirect 
exchange coupling between two semi-infinite ferromagnetic strips (FM electrodes) separated by 
metallic graphene nanoribbons (GNRs) is investigated. In order to form a FM/GNR/FM junction, 
>— ^ ' a graphitic region of finite length is coupled to the FM electrodes along graphitic zigzag or armchair 

interfaces of width N. The numerical results show that, the exchange coupling strength which can 
^ . be obtained from the difference between the total energies of electrons in the ferromagnetic and 

antiferromagnetic couplings, has an oscillatory behavior, and depends on the Fermi energy and the 
| length of the central region. 

vo : 

^ ! I. INTRODUCTION 

Magnetic exchange interaction between ferromagnetic (FM) layers separated by a nonmagnetic metal spacer has 
been investigated to a large extent because of both possible applications exploiting the giant magnetoresistance effect 
Cj and the expectation of the way to the production of nonvolatile computer memories, extremely efficient magnetic 
sensors, and magnetic materials with enhanced information storage capacity 0, Among low-dimensional struc- 
tures, graphene nanoribbons (GNRs) and carbon nanotubes (CNTs) are potentially useful for magneto-transport 
applications. For that purpose, GNRs and CNTs must interact with magnetic foreign objects such as substrates, 
impurities, adsorbed atoms, and nanoparticles [!, BJ. When the separation between the magnetic objects increases, 
'■^ ■ the direct exchange interaction between magnetic moments decays abruptly since this interaction requires a finite 
overlap between the wave functions that surround the respective magnetic objects. Therefore, when the moments 
are not too close together there is no direct overlap between their wave functions. In such a case, the only way a 
i magnetic coupling can arise is if it is mediated by the conduction electrons of the metallic host. Understanding the 

physics of this so-called indirect exchange coupling between localized magnetic moments mediated by the conduction 
electrons of metallic hosts might provide clues to making magnetic carbon-based structures which could be applicable 
^ | in spintronic devices 

In the field of exchange coupling of magnetic trilayers or multilayers, many methods and models have been employed, 
such as the first-principles method @, the tight-binding total energy calculation 0], the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) theory the one-band tight-binding hole-confinement model {9\, the free-electron model etc. 
In this article, based on the single-band tight-binding approximation and the self-energy method, we investigate the 
indirect exchange coupling between two semi-infinite FM electrodes separated by metallic armchair GNRs (AGNRs) 
qq ■ and zigzag GNRs (ZGNRs). 
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II. MODEL AND FORMULATION 

x. 

We consider a magnetic nanostructures consisting of two semi-infinite FM nanostripes (electrodes) separated by a 
graphitic region of length L (in the x direction) and width N (in the y direction). Since we study the exchange coupling 
between two FM electrodes, the electronic structure of the central part of the junction (i.e. GNR) should be resolved 
in detail. It is therefore reasonable to decompose the total Hamiltonian of the system as H = Hl + Hr + He + V. For 
simplicity, we treat a square lattice structure for the FM electrodes. The Hamiltonian for such electrodes is written 
within the tight-binding approximation as 

= 51 ' h «)C<T^»,<T - *Wc,C,<A*, CT , (1) 
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where cj creates (destroys) an electron with spin er at site j in electrode a (=L, R), and U a j a = t is the 

hopping matrix element between nearest-neighboring sites « Q and j a . Here, e Q is the spin independent on-site energy 
and will be set to zero, — a ■ h a is the internal exchange energy with h a denoting the molecular field at site i a , and 
a being the conventional Pauli spin operator. On the other hand, the Hamiltonian of GNR in the absence of FM 
electrodes is expressed as 



E 

<1C,3C> 



(2) 



where d\ c a (di Cy(T ) creates (destroys) an electron with spin a at site i of GNR. Finally, V describes the coupling 

between the electrodes and the GNR and takes the form V = — J2i a j c U a ,j c (^l a a^io," + H.c). The hopping 
elements ti c Jc between the tt orbitals of the GNR, and also t ia j G between the lead orbitals and the tt orbitals of the 
GNR are taken to be i; that is, the hopping parameters are same in the three regions. 

In this study we assume that the spin direction of the electron is conserved in the propagating process through 
the GNR. Therefore, there is no spin-flip process and the spin-dependent transport can be decoupled into two spin 
channels: one for spin-up and the other for spin-down. On the other hand, it is well known that the effect of semi- 
infinite electrode a on the GNR can be described by a (spin-dependent) self-energy matrix £ a ,<r llj. Therefore, it is 
reasonable to write H = '}2 IJ H r j, where H a = He + £l.<t + £r,ct- Now the spin-dependent Green's function of the 
GNR coupled to the two FM electrodes is given as 



GW(e) = [el -H c - E i , (T (e) - S fl , CT (e)] 



(3) 



where the self-energy matrices contain the information of the electronic structure of the FM electrodes and their 
coupling to the GNR. These matrices can be expressed as 



S a ,o-(e) = Tc, a gu,a{t)T a ,c 



(4) 
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FIG. 1: J as a function of AGNR length for E F =0.5 
eV. 



FIG. 2: J as a function of ZGNR length for E F =0.5 
eV. 



where f is the hopping matrix that couples the GNR to the leads and is determined by the geometry of the GNR- 
lead bond. g a ^ are the surface Green's functions of the uncoupled leads i.e., the left and right semi- infinite magnetic 
electrodes, and their matrix elements are given by 

4 ^ sm(fca;i)sin( ] ^y i )sin(A;a; i )sin( ] ^ T j/j) 
g a Ahj;z) - Nx{Ny - 1} 2_, - _ €a - - . ha _ 2t[cos{ka) + 008(^)1 (5) 

Here, z — e + iS, 1 < I < N y , k £ [— f , f ] and Np with /3 = x, y is the number of lattice sites in the j3 direction. In 
order to calculate the total energy, we need the total density of states per site for the electron of spin a (=f or j) 
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which is given as D a (e) = —-^ Yli I m [G<j(e)]u where 3\f is the total number of carbon atoms in the GNR. Thus, at 
T = K, the total energy of the electrons that occupy the levels up to the Fermi energy ep in FM (antiferromagnetic 
(AF)) configuration is given by 



£™W = f e[D FM{AF \e) + D FM{AF \e)]de 

J — OG 



(6) 



The indirect exchange coupling J, is defined as the total energy difference between the FM and AF configurations 
of the system: 



j = gFAf _ £ AF 



(7) 



Based on the above formalism, we calculate J for two different junctions 
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FIG. 3: J as a function of AGNR length for E F =1 
eV and E F =2 eV. 
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FIG. 4: J as a function of ZGNR length for E F = 1 
eV and E F =2 eV. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



We have done the numerical calculations for the case that the direction of magnetization in the left FM electrode is 
fixed in the +y, while the magnetization in the right electrode is free to be flipped into either the +y or -y direction. 
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We set |h a |=1.5 eV and t=l cV in the calculations. Our results clearly show that J oscillates with increasing the 
length of graphitic region in both FM/AGNR/FM (Figs.l and 3) and FM/ZGNR/FM (Figs.2 and 4) junctions. It 
can be seen from the figures that as the length of GNR increases, the oscillation amplitudes of J decrease. Also, 
the oscillatory phase changes by varying the Fermi energy. We can ascribe the indirect coupling to the interferences 
of electron waves in the metallic GNRs because of spin-dependent reflections at the FM/GNR interfaces. Electron 
states in the GNR can be propagative Bloch waves, but also evanescent waves due to the finite length of the spacer. 
Propagative states give rise to oscillatory contributions to the coupling, while evanescent states yield damping terms. 

Therefore, both kinds of states contribute to J. Their importance depends essentially on the nature of the states at 
the Fermi level. It is clear from the figures that, in most cases, when the separation between two electrodes exceeds 
15 unit cells, the coupling between two magnetic electrodes becomes very small. These results may be useful for 
understanding the strength of exchange coupling between two magnetic impurities or adsorbed atoms on GNRs and 
CNTs. 
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